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C 2 ^ ' Geometric phases accompanying adiabatic processes in quantum systems can be 

, utilized as unitary gates for quantum computation. Optimization of control of the 

C ~_ ~) • adiabatic process naturally leads to the isoholonomic problem. The isoholonomic 

I/"") ' problem in a homogeneous fiber bundle is formulated and solved completely. [Pro- 

C " 3 ' ceedings of International Conference on Topology in Ordered Phases organized by 

, Hokkaido University in March 2005.] 

-i— > ■ 

£^ ' 1- Introduction 

& '. 

^ i The isoholonomic problem was proposed in 1991 by a mathematician, 

CT 1 - Mont gomery . The isoholonomic problem is a generalization of the isoperi- 

metric problem, which requests finding a loop in a plane that surrounds the 
largest area with a fixed perimeter. On the other hand, the isoholonomic 
problem requests finding the shortest loop in a manifold that realizes a 
specified holonomy. This kind of problem naturally arose in studies of the 
Berry phase 2-4 and the Wilczek-Zee holonomy 5 , which appear in a state of 
a controlled quantum system when the control parameter is adiabatically 
changed and returned to the initial value. Experimenters tried to design 
efficient experiments for producing these kinds of holonomy. Montgomery 
formulated the isoholonomic problem in terms of differential geometry and 
gauge theory. Although he gave partial answers, construction of a concrete 
solution has remained an open problem. 

Recently, in particular after the discovery of factorization algorithm by 
Shor 6 in 1994, quantum computation grows into an active research area. 
Many people have proposed various algorithms of quantum computation 
and various methods for their physical implementation. Zanardi, Rasetti 7 
and Pachos 8 proposed utilizing the Wilczek-Zee holonomy for implementing 
unitary gates and they named the method holonomic quantum computa- 
tion. Since holonomy has its origin in geometry, it dose not depend on 



1 
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detail of dynamics and hence it does not require fine temporal tuning of 
control parameters. 

It should be noted, however, that holonomic quantum computation re- 
quires two seemingly contradicting conditions. The first one is the adi- 
abaticity condition. To suppress undesirable transition between different 
energy levels we need to change the control parameter quasi-stationarily. 
Hence a safer control demands longer execution time to satisfy adiabaticity. 
The second one is the decoherence problem. When a quantum system is 
exposed to interaction with environment for a long time, the system loses 
coherence and a unitary operator fails to describe time-evolution of the 
system. Hence a safer control demands shorter execution time to avoid de- 
coherence. To satisfy these two contradicting conditions we need to make 
the loop in the control parameter manifold as short as possible while keep- 
ing the specified holonomy. Thus, we are naturally led to the isoholonomic 
problem. 

We would like to emphasize that a quantum computer is actually not 
a digital computer but an analog computer in its nature. Hence, the ge- 
ometric and topological approaches are useful for building and optimizing 
quantum computers. 

This paper is based on collaboration with D. Hayashi and M. Nakahara 9 . 
We are further developing our studies on optimal and precise control of 
quantum computers with Y. Kondo, K. Hata and J.J. Vartiainen 10-12 . We 
thank Akio Hosoya, Tohru Morimoto and Richard Montgomery for their 
kind interest in our work. 

2. Wilczek-Zee holonomy 

A state vector ip(t) £ C N evolves according to the Schrodinger equation 

mj t m = mm- « 

The Hamiltonian admits a spectral decomposition H(t) — ^2f = i£i(t)Pi(t) 
with projection operators Pi(t). Therefore, the set of energy eigenval- 
ues (ei, . . . , el) and orthogonal projectors (Pi, . . . , Pi) constitutes a com- 
plete set of control parameters of the system. Now we concentrate on the 
eigenspace associated with the lowest energy E\. We write Pi(t) as P{t) for 
simplicity. Suppose that the degree of degeneracy k = tr P(t) is constant. 
For each t, we have the eigenvectors such that 



H(t)v a (t)=e 1 (t)v a {t), (a = !,...,*;). (2) 
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We assume that they are normalized as vl,(t)vp(t) = 5 a p. Then 

V(t)= («i (*),..., « fc (t)) (3) 

forms an N x k matrix satisfying V^(t)V(t) = I k and V{t)V\t) = P(t). 
Here I k is the fc-dimensional unit matrix. The adiabatic theorem guarantees 
that the state remains the eigenstate associated with the eigenvalue £\(t) 
of the instantaneous Hamiltonian H (t) if the initial state was an eigenstate 
with £i(0). Therefore the state vector is a linear combination 

k 

= £><,(*) «a(t) = v(*M*)- (4) 

a=l 

The vector <p = t {(j>\, . . . ,(f>k) G C fe is called a reduced state vector. By 
substituting it into the Schrodingcr equation (1) we get 

Its solution is formally written as 

<Kt)=exp(-±J ei (8)(k)Texp(- J yt ^r ds )^°)' ( 6 ) 
where T stands the time-ordered product. Then ip(t) = V(t)cj>(t) becomes 
0(t) = e -Ho«(«)^y( t )Te-^ ytdy y t (0)V(0). (7) 

In particular, when the control parameter comes back to the initial point as 
P(T) = P(0), the state vector ip(T) also comes back in the same eigenspace 
as -0(0) = V (0)0(0). The Wilczck-Zcc holonomy T G U(k) is defined via 

i>(T) ^e-^ lo £l ^ ds V(0)T(t>(0) (8) 

and is given explicitly as 

r = V(0)*V(T)Te-f v * dv . (9) 

If the condition — is satisfied, the curve V{t) is called a horizontal 

lift of the curve P(t). Then the holonomy (9) is reduced to T = V^(0)V(T). 

3. Formulation of the problem 

The isoholonomic problem is formulated in terms of the homogeneous fiber 
bundle (SW,fc(C), G Nik (C),ir,U(k)). The Sticfcl manifold S Ntk (C) is the 
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set of ortho-normal fc-frames; a fc-frame V spans the degenerate energy 
eigenspace in C^; 

S N ,k(Q = {V G M{N, k; C) | V*V = h}, (10) 

where M(N, k; C) is the set of N x k complex matrices. An element of the 
unitary group h 6 U(k) acts onFe Sjv,fc(C) from the right as (V, /i) Vh 
by means of a matrix product. The Grassmann manifold GN,k(C) is defined 
as the set of projection matrices to fc-dimensional subspaces in C N , 

G N , k (C) = {P e M(N,N;C) \ P 2 = P, P f = P, trP = k}. (11) 

The projection map tt : 5jy.fc(C) — > Gjv,fc(C) is defined as 7r : K P := 
FT^. Then it can be proved that the Stiefel manifold <SW,fc(C) becomes a 
principal bundle over Gjy j fc(C) with the structure group ?7(fc). The canon- 
ical connection form on 5iv,fc(C) is defined as a one- form 

A = VUV, (12) 

which takes its value in the Lie algebra u(fc). The holonomy associated 
with this connection is called the Berry phase in case of k = 1 and 
the Wilczek-Zee holonomy in general. We define Riemannian metrices, 
\\dV\\ 2 = tr(dVUV) for the Stiefel manifold and ||dP|| 2 = tr (dPdP) for 
the Grassmann manifold. For any curve P(t) in Gjv,fe(C), there is a curve 
V(t) in S N . k (C) such that 7r(V(i)) = P(t). If the curve V(t) satisfies 

(.3) 

it is called a horizontal lift of the curve P(t). When the curve P(t) is a 
closed loop, such that 

V{T)V\T) = V(0)V\0), (14) 

the holonomy associated with the loop is defined as V(T) = V(0)T and is 
given as 

r = vH0)V(T) e U(k). (15) 

We formulate the isoholonomic problem as a variational problem. The 
length of the horizontal curve V(t) is evaluated by the functional 

where ft(t) E u(k) is a Lagrange multiplier to impose the horizontal con- 
dition (13) on the curve V(t). Thus the isoholonomic problem is stated as 
follows; find a horizontal curve V(t) that attains an extremal value of the 
functional (16) and satisfies the boundary conditions (14) and (15). 
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4. Derivation and solution of the Euler-Lagrange equation 

We derive the Euler-Lagrange equation associated the functional S and 
solve it explicitly. A variation of the curve V(t) is defined by an arbitrary 
smooth function r](t) £ u(N) such that r/(0) = r](T) = and an infinitesimal 
parameter e £ R as 

V e (t) = (l + e V (t))V(t). (17) 

By substituting V e (t) into (16) and differentiating with respect to e, the 
extremal condition yields 

J trj^VVt -VW -VQV^dt. (18) 



0= ^ 

de 



e=0 

Thus we obtain the Euler-Lagrange equation 

j(VV* -VV ] + VW t ) = 0. (19) 

The extremal condition with respect to Ct(t) reproduces the horizontal equa- 
tion V*V = 0. 

Next, we solve the equations (13) and (19). The equation (19) is inte- 
grated to yield 

VV* -VV^ + VSlVi = const = X £u(N). (20) 

Conjugation of the horizontal condition (13) yields V^V = 0. Then, by 
multiplying V on (20) from the right we obtain 

v + vn = xv. (21) 

By multiplying V* on (21) from the left we obtain 

Q = V^XV. (22) 

We can show tl = by a straightforward calculation. Hence, O(t) is actually 
a constant matrix. The solution of (21) and (22) is 

V(t) = e tx V e~ ta , n = VjXV . (23) 

We call this solution the horizontal extremal curve. Then (20) becomes 

(xv - vnjvt - v(-v^x + nv r ) + vnv 1 = x, 

which is arranged as 

X - (W^X + XW ] - W ] XVV^) = 0. (24) 
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Here we used (22). We may take, without loss of generality, 

Vo=( 1 ,':) eS N .k(C) (25) 



^ , 

as the initial point. We can parametrize X <G u(N), which satisfies (22), as 

with W G M(k, N-k;C) and Z G u(N - k). Then the constraint equation 
(24) implies that Z = 0. Finally, we obtained a complete set of solution 
(23) of the horizontal extremal equation (13) and (19). 

5. Solution to the boundary value problem 

The remaining problem is to find the controller matrices £1 and W that 
satisfy the closed loop condition 

V{T)V\T) = e TX V V^e- TX = V V^ (27) 

and the holonomy condition 

V(T) = e TX V e~ Tn = C/ gatc (28) 

for a requested unitary gate U gate G U(k). Montgomery 1 presented this 
boundary value problem as an open problem. Here we give a prescription 
to construct a controller matrix X that produces the specified unitary gate 
f/gate- It turns out that the working space should have a dimension N > 2k 
to apply our method. In the following we assume that N — 2k. The time 
interval is normalized as T = 1. 

Our method consists of three steps. In the first step, we diagonalize a 
given unitary matrix f7g atc G U(k) as 

B)U &ate R = U diag = diag(e 4 ^ , . . . , e 4 ^ ) (0 < j 3 < 2ir) (29) 

with R G U(k). The small circle is a circle in a two-sphere CP 1 C Gjv,fe(C) 
that surrounds a solid angle which is equal to twice of the Berry phase. In 
the second step, combining k small circles we construct k x k matrices 

O diag =diag(iwi,...,iw fe ), W d i ag = diag(in, . . . ,ir k ), (30) 

with ujj = 2(7r — 7j) and Tj = yV 2 — (n — 7j) 2 - We combine them into 
a, 2k x 2k matrix 



X d 



lag 



^diag ^diag 
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In the third step, we construct the controller X as 

f^diag TT diag ^ ftf \ _ ( i?r! diag i?t RW diag 



X 





Oh) \-WlR 1 



(31) 

In the paper 9 we calculated explicitly controllers of various unitary gates; 
the controlled NOT gate, the discrete Fourier transformation gate and so 
on. 



6. Conclusion 

We formulated and solved the isoholonomic problem in the homogeneous 
fiber bundle. The problem was reduced to a boundary value problem of the 
horizontal extremal equation. We determined the control parameters that 
satisfy the boundary conditions. This result is applicable for producing 
arbitrary unitary gates. 
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